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Abstract 

The Turaev-Viro invariant for a closed 3-manifold is defined as tfie contraction of a certain 
tensor network. The tensors correspond to tetrahedra in a triangulation of the manifold, with 
values determined by a fixed spherical category. For a manifold with boundary, the tensor 
network has free indices that can be associated to qudits, and its contraction gives the coefficients 
of a quantum error-correcting code. The code has local stabilizers determined by Levin and Wen. 
For example, applied to the genus-one handlebody using the Z2 category, this construction yields 
the well-known toric code. 

For other categories, such as the Fibonacci category, the construction realizes a non-abelian 
anyon model over a discrete lattice. By studying braid group representations acting on equiva- 
lence classes of colored ribbon graphs embedded in a punctured sphere, we identify the anyons, 
and give a simple recipe for mapping fusion basis states of the doubled category to ribbon graphs. 
We explain how suitable initial states can be prepared efficiently, how to implement braids, by 
successively changing the triangulation using a fixed five-qudit local unitary gate, and how to 
measure the topological charge. Combined with known universality results for anyonic systems, 
this provides a large family of schemes for quantum computation based on local deformations 
of stabilizer codes. These schemes may serve as a starting point for developing fault-tolerance 
schemes using continuous stabilizer measurements and active error-correction. 
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1 Introduction 



Proposed topological quantum computers protect quantum information in a physical medium, 
while allowing logical operations to be applied robustly using adiabatic local processes [Kit03, 
Pre04, FKLW03, NSS"'"08]. Building systems with the required properties is generally believed 
to be difficult. However, some of the highest fault-tolerance threshold estimates are based on a 
code — Kitaev's toric code — and procedures that are motivated by topological quantum computa- 
tion [RHG06, RH07, DKLP02, FSG09]. A limitation of the toric code is that it corresponds to an 
abelian anyon model, so does not allow for universal quantum computation. By translating more 
general anyon models into the language of error-correcting codes, this limitation can be surpassed. 

A first step is to realize the Hilbert space of a topological quantum field theory (TQFT) on 
a spin lattice. Turaev and Viro completed this for a class of ribbon categories [TV92], and their 
construction was extended by Barrett and Westbury to spherical categories [BW99, BW96]. For 
example, if C is the Fibonacci category, a modular tensor category with one non-trivial particle r, 
the Turaev- Viro construction gives a TQFT for the "doubled" category C <^ C* , with particles 
l(S'l,l(8)r, riSil and t^t. Levin and Wen showed that the TQFT Hilbert space is the code space 
of a set of commuting local stabilizers [LW05b]. 

Roughly, a TQFT associates to every surface S a Hilbert space 1-Ly,-, and to every diffeomorphism 
/i : S — )■ S' a linear map U{h) : T-Ly. — ^ ^s'- By restricting to self-diffeomorphisms of S, the maps 
{U{h)}h form a representation of S's mapping class group. For example, if S is the sphere with 
n punctures, its mapping class group is the braid group on n — 1 strands. For certain TQFTs, 
including that based on the Fibonacci category, the resulting computational model has been shown 
to be equivalent to standard quantum computers [FLW02, FKW02]. 

We address the problem of implementing TQFT computations for modular tensor categories in 
the Turaev- Viro spin-lattice code space. This involves the following three steps: 

(i) We identify bases of T-Ly. corresponding to fusion diagrams for the doubled theory. This 
essentially boils down to decomposing the braid group representation on T-Ly,^ a Hilbert space 
of colored ribbon graphs modulo local equivalence relations. 

(ii) For every braid group element 6, we show how to realize the encoded unitary U{b) as a sequence 
of protected, local gates. Pachner moves on the surface's triangulation correspond to unitary 
F- moves on the spins. Under such a transformation, the code space, i.e., the ground space Ht 
of a local stabilizer Hamiltonian Ht, is transformed to the ground space of a local stabilizer 
Hamiltonian Ht+i which can be obtained by adiabatically changing three terms in Ht- 

(iii) Finally, we give procedures for preparing certain initial states in T-Ly^ and for measuring the 
topological charge within a region. 

Altogether, this provides a large family of schemes for universal quantum computation using geo- 
metrically local operations on locally stabilized codewords. These schemes may serve as a starting 
point for developing fault-tolerance schemes using continuous stabilizer measurements and active 
error-correction. 

The computational scheme presented here is intimately linked to fundamental concepts in topo- 
logical quantum computation. The existence of such a computational scheme based on a doubled 
category was previously conjectured in [RSW09, Section 6.2]. A similar scheme was proposed by 
Freedman using the SU{2) Witten-Chern-Simons modular functor at level three [FreOO], but in 
contrast to our work the stability of his proposal relies on a conjectured energy gap. 
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The purpose of this paper is partly introductory and we keep the presentation as self-contained 
as possible. In particular, though helpful, we will not require previous knowledge of anyonic fusion 
spaces; see [Pre04] and the appendix of [Kit06] for excellent introductory reviews. 

The ideas in this paper originated in an email exchange between one of the authors (G.K.) and 
Alexei Kitaev. 



Outline 

In Section 2, we introduce the Hilbert space "Hs associated with a surface E by the Turaev-Viro- 
TQFT defined by a category C. This is a Hilbert space of colored ribbon graphs embedded in S 
modulo local equivalence relations. We explain how self-diffeomorphisms /i : S — )• S act on this 
space by deforming these graphs. In Section 3, we show how to decompose this mapping class group 
action into its anyonic content: the representation is described by the "doubled" category C (dC*. 
In Section 4, we explain how the space T-Ls arises as a subspace of qudits on a lattice; this is the 
Turaev-Viro code. We also introduce the corresponding stabilizer Hamiltonian defined by Levin 
and Wen [LW05b]. In Section 5, we show how to encode information into the code space, execute 
logical operations, and read out encoded information. 

Throughout Sections 2 through 5, we restrict our attention to the Fibonacci category. In 
Section 6, we extend these results to arbitrary ribbon categories. We conclude in Section 7 with a 
discussion of the relation between the Turaev-Viro 3-manifold invariant and the Turaev-Viro code. 

In the appendices, we give proofs and examples for various claims, and show how our results 
extend to mapping class group representations of surfaces with higher genus. In particular, this leads 
to an alternative (known) expression for the Turaev-Viro invariant, as discussed in Appendix C. 



2 The Hilbert space I-Ly. of Fibonacci ribbon graphs on E 
2.1 Definition of Fibonacci ribbon graphs and He 

Let S be a compact, orientable surface with boundary, and a single marked point on each boundary 
component. We begin by associating to S a "ribbon graph" Hilbert space ^s- 

A ribbon graph (also called skein/spin network or string net) on S is the smooth embedding 
into the interior of S of a graph in which each vertex has degree either two or three, except that a 
single vertex of degree one is allowed to be mapped to each marked point on a boundary component. 
The ribbon graph Hilbert space T-Ly; is the space of formal linear combinations of ribbon graphs 
embedded in S modulo the following local relations: 




Here the first rule is meant to indicate that ribbon graphs related by smooth deformations of the 



4 



embedding are equivalent. The second rule allows us to eliminate loops, picking up a factor of 

r:=^, (4) 

provided that the loop can be contracted to a point within the surface. The third rule allows for 
eliminating adjacent trivalent vertices from the ribbon graph. 
By applying the above rules, we derive 

M3 ) ) O =0 (5) 

Therefore we impose as an additional rule 



^ =0 (6) 

The only case in which this is not already a consequence of the other rules is when the dangling 
edge connects to a boundary point. 



2.2 Computational bases for and inner product on 

The space T-L-£ is finite dimensional, with a dimension that depends on the topology of S. We will 
concentrate on the case where S is the n-punctured sphere, S„ = S*^ \ (A^ U • • • U A^), i.e., the 
sphere with n discs A^ removed and one boundary point pi G dA^ for every hole. Then 
has some simple bases associated to the dual graphs of certain triangulations. We call these bases 
"computational bases" because they are useful for encoding 7iY;„ using qubits. 

A labeling ^ of S associates an element £{p) G {0, 1} to every marked boundary point p. A 
ribbon graph is consistent with i if it has an edge ending at exactly those p with i{p) = 1. Then 
the space "Hs can be decomposed as 

^E = 0^i, (7) 

e 

where is the subspace spanned by ribbon graphs that are consistent with the labeling £. 

First consider the case where we restrict the ribbon graph to have open boundary conditions, 
i.e., to have no edges touching the boundaries of S. In this case, fill in every boundary piece with 
a disk, and consider a (degenerate) triangulation of S whose vertices are the centers of these disks. 
Then the graph dual to this triangulation is trivalent. Placing a bit, or 1, in the center of each 
edge, determines a ribbon graph by interpreting a 1 as the presence of an edge and a as an absence. 
The set of all such 0/1 assignments that have no dead ends — that is, for every vertex of the dual 
graph, zero, two or three of the incident edges must be set to 1 — is a basis for "He- Moreover, we 
can define an inner product on Tis by setting these basis states to be orthonormal. One can verify 
that this inner product does not depend on the choice of basis, i.e., on the triangulation. Several 
examples are given in Figure 1. 

The fully general setting, in which ribbon graph edges are allowed to go to a marked point on 
each boundary component, is slightly more complicated. In this case, again fill in the holes and 
place a vertex in the middle of each hole. For each of these vertices, add a loop going around the 
marked point on the boundary. Then complete this graph to a triangulation of the vertices. The 
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(a) (b) (c) (d) 



Figure 1: In (a) and (b) are shown two different examples of basis choices for T-L^, for the case of S 
being the sphere with three punctures (indicated with crosses). Dashed hues mark the point-set 
triangulation of the punctures, while the dual graph, which is trivalent, is indicated with solid lines. 
Part (c) gives another example of a basis for the sphere with four punctures. In (a-c), it is assumed 
that there are no marked points on the boundaries of S. Part (d) shows the more general situation, 
for the sphere with two shaded holes, each of which has a marked point on its boundary. 

dual graph will be trivalent, except with dead ends at the marked boundary points. Then the basis 
and inner product are defined as above on this graph. 

For example, if n = or 1, then T-Ly.„ = C is one-dimensional. (In the case n = 1, the 
rule in Eq. (6) ensures that T-L^^, ribbon graphs with an edge attached to the boundary, is zero- 
dimensional.) Figure 1(d) shows the n = 2 case. Note that there are four edges in the dual 
graph. Since each edge can be absent or present (0 or 1), this gives a natural isometry of Ht.a 
into (C^)®^. However, there are only seven edge assignments without dead ends; dim(^s4) = 
dim(-?^|';°^)+dim(^g;^^)+dim(-H|';^^)+dim(1^^^;°^) = 2 + 3+1 + 1 = 7. More generally, for n > 2, 
a basis derived in this way from a point-set triangulation gives a natural isometry from T-Ly.„ into 
^^2^®{5n-6)^ i.e., a way of implementing T-Lt.^ using 5n — 6 qubits. Therefore, we call such a basis 
for T-Ly.^ a computational basis. It defines an inner product on • 

The examples in Figure 1 suggest a general procedure for defining a basis for , based on a 
pants decomposition of the surface S„. For example, the pants decomposition of S4 corresponding 
to Figure 1(c) is 




(8) 



In general, the pants decomposition corresponds to a rooted binary tree with n — 1 leaves, the root 
and each leaf associated to one of the holes in The degree-three internal vertices correspond to 
"pants segments," each isomorphic to S2, and the edges are "cylindrical segments," each isomorphic 
to S3. The computational basis is obtained by associating four qubits to every leaf, and one qubit 
to the root and to every internal pant segment. 
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2.3 Actions on "Hs 



Let and be two surfaces, and let S be the surface obtained by gluing together one or more 
boundary circles of the surfaces, in such a way that the marked boundary points are matched. 
Then ribbon graphs on Tia and T,b with matching boundary labels can be glued together to obtain 
a ribbon graph on S. In this way, every element s G defines a linear operator s : 'H-£j^ ^s- 
The result of applying s to t £ T~Lt.a the linear combination of ribbon graphs obtained by gluing 
together all ribbon graphs comprising s and t with matching labels, up to the equivalences in 
Section 2.1. The map ^ : T-Lt,b ~^ clearly an isomorphism of vector spaces; we will therefore 

interpret ribbon graphs interchangeably as states or operators. For example, if = S2 and 
= S3, 




In alternative notation, the same equation can be written 

9 (« • Y + ^ ■ ) = • 

This result can be expressed in terms of the basis ribbon graphs by application of the rules (l)-(3). 

Now specialize once more to the case of the n-punctured sphere S = The surface S together 
with its marked boundary points pi, ■ ■ ■ ,Pn defines the mapping class group A^(S, {pi, ■ ■ ■ ,Pn})- Its 
elements are isotopy classes of orientation-preserving diffeomorphisms of S which fix the boundary 
points. The mapping class group is generated by operators of two kinds, Dehn-twists and braid- 
moves. A Dehn-twist is a 27r-twist along a simple closed curve. A braid-move affects a pair of pants 
with two holes in S, and is defined as a vr- twist around a simple closed curved on S enclosing both 
holes, followed by (or preceded by) vr-twists on each of the legs, viz. 




The mapping class group Ai acts on the ribbon graph space by linearly extending its action 
on the basis ribbon graphs to the full space Tij:. This action is unitary with respect to the inner 
product defined previously. 

Let us give a more explicit description of the action of Dehn-twists and braid-moves on the 
ribbon graph Hilbert space. First consider a non-contractible curve 7. 




(12) 



The surface S is topologically unchanged if we cut along 7 and insert into the gap a cylindrical 
segment S2. Then the Dehn-twist about the 7 corresponds to the linear operator .0(7) : T^s — ^ 
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given by 

= ? + 9 (13) 



Similarly, the operator corresponding to braiding can be specified by 



B = 




where the binary variables i and j indicate the presence or absence of an edge. Indeed, stacking 
this on top of the two legs in the diagram on the left-hand side of Eq. (11) gives the right-hand 
side thereof. 

We will next identify irreducible subspaces of the mapping class group action. 




3 Doubled Fibonacci anyons in 

Let S„ be the n-punctured sphere. Section 2.2 showed how to construct a computational basis for 
^s„ based on one of a family of point-set triangulations. It also specialized this construction to the 
case of a triangulation derived recursively from a pants decomposition of S„. These computational 
bases are useful for encoding using qubits, and we will study this implementation further in 
Section 4. 

However, these bases are inconvenient operationally, because they do not transform cleanly 
under the application of Dehn-twists and braid-moves. In particular, it is complicated to express 
in the pants decomposition computational basis the action of a Dehn-twist about either the root 
or an internal cylindrical segment. From the right-hand side of Eq. (10), the loop has to be 
propagated upward toward the leaves, potentially requiring many applications of the ribbon graph 
equivalence rules. Additionally, the basis is asymmetrical, since the puncture at the root of the 
pants decomposition tree is treated differently from the other punctures. Intuitively, too much 
information is kept at the leaves. 

In this section, we will define an orthonormal basis that simultaneously diagonalizes the Dehn- 
twists. Each basis element will be represented by an "anyon fusion diagram." Although slightly more 
complicated to define, the advantage of this anyonic basis is that it will reveal more mathematical 
structure and greatly clarify the steps necessary for obtaining universal quantum computation. 

Let a bold, dashed line indicate a linear combination of the two diagrams with and without 
that line: 

^ ( +r ) (15) 



Such lines are called vacuum lines, because a simple calculation shows that other lines can pass 
over them freely: 




vi+^\(0) = \ 0+A(0)= (m = Vi + T^ (m (16) 
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Here, Eq. (16) is to be interpreted as a local identity, that is ribbon graphs inside the hatched area 
or outside the diagram are assumed to be the same for every term. Its proof is a straightforward 
application of the rules (l)-(3). Lemma A.l states further properties of the vacuum lines. 



3.1 Anyonic fusion basis for "Hsn 

Fix a pants decomposition of S„, corresponding to a rooted binary tree T. Fix a labeling £ of the 
marked boundary points, as in Eq. (7). 

A labeling of T is an assignment to each edge of T either 1 or r. A labeling is fusion-consistent 
if no internal vertex has exactly one incident edge labeled r. A pair of labelings of T is boundary- 
consistent with i if for each boundary point p, with labels of the corresponding edges 5+ and 
b-, 

£(p) = 0^6+6_e {11, tt} 

i{p) = l^b+b-e{Tl,lr,rT} . 

A pair of fusion-consistent labelings of T which is boundary-consistent with i is called an ^-consistent 
doubled anyon fusion diagram. 

The anyonic fusion basis states for Ti^^^ are an orthonormal basis, with basis elements indexed 
by ^-consistent doubled anyon fusion diagrams. In the remainder of this subsection, we will define 
the state \i,d) G that corresponds to an ^-consistent doubled anyon fusion diagram d. The 

state \£, d) is defined by first constructing a three-dimensional ribbon graph living in the thickened 
surface S x [—1,1], and then reducing this ribbon graph down to two dimensions. 

Think of each of the two labelings in d as a ribbon graph, by interpreting label 1 as the absence 
of an edge and label t as the presence of an edge along the edges of the tree T. Place these two 
ribbon graphs on S x {1} and Sx{— 1} inSx [—1,1]. Then add vacuum lines around each puncture 
in E X {0}, except not around the puncture at the root of the pants decomposition. 

In the resulting three-dimensional ribbon graph, edges can end at (p, — 1) and (p, 1) for each 
marked boundary point p of S. However, we would like there only to be a single boundary condition 
at the point (p, 0). To fix the boundary conditions, close off the ribbon graphs at each boundary 
with one of the five following diagrams. Here, the vertical axis represents the second coordinate, 
from —1 to 1, with the point (p, 0) marked. 




Boundary-consistency of d implies that one of these diagrams applies. Note that this closing off is 
done inside the vacuum line around the hole at p. 

Finally, visualize the three-dimensional ribbon graph as a two-dimensional diagram with cross- 
ings. (This requires slightly offsetting the ribbon graph boundaries in S x {1} from the boundaries 
in E X {—1}. Follow any convention for choosing the direction of the offset; different conventions 
will only change the phase of \i, d).) Use repeatedly the rule 

\/ ^ g-3W5 ) ( + e3-/5 ^ (19) 
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to eliminate ribbon crossings, finally obtaining a ribbon graph in S x {0} = S. This state is \i, d). 

Note that ribbon graphs in three dimensions can be manipulated similarly to ribbon graphs 
in two dimensions, for example by smoothly deforming the embedding. However, it is useful to 
think of the ribbon graphs in three dimensions as actually being thin ribbons, since kinks cannot 
be freely eliminated. For example, we derive from Eq. (19) that 



-A-Ki 



/5U 



(20) 



The anyonic fusion basis states form an orthonormal basis. We will prove this and other 
properties, in a more general setting, in Appendix A. 2. 

Example: Anyonic fusion basis states for 

For T,2, the only pants decomposition is as a single cylindrical segment. T-Lj:^ is seven-dimensional, 
spanned by the states {;, 1,9,9,0,0,9}, defined by 



(21) 
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o 
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To explain the notation, notice, for example, that 9 corresponds to looking downward on the final 
cylinder above, as the base of the cylinder spreads out. The basis elements are all of unit length 
and orthogonal, except (9I9) = 1/t. 

An anyonic fusion basis for Ti^^ is determined by a single label 11, tt, It or rl consistent 



with the boundary conditions £. Explicitly, the orthonormal bases for Hj^^ = ^ 
^,(1,0) ^ ^,(0,1) . , 



(0,0) 

S2 



n 



(1,1) 
S2 



n 

n 
n 



=(0,0) 



^=(1,1) 



^=(1,0) 
^=(0,1) 

S2 




vr+7211 



+ rro 

5^/'--i|,/, 
\/1 + t2 It 



a/1 + T2rl 

{\£,tt) :=|9)} 
{\e,TT) :=|9)} 




(22) 



Here 



T : = 



1 

T 



1 - Vb 



(23) 
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and we are using the five operators 



^^o = rT?2('^+^°) -^^1 = 7^(^7 + 9 + 9) (24) 

1 1 

It = Y^~^ ( 'i^ + ^^"^ 9 + ™^ 9 ) T"l = Yip72 ( T + ™^ 9 + rco'^ ^ ) 

with w = e'^*/^. These operators are the five minimal central idempotents in the algebra of operators 
T~{-T:2 ~^ of the form s, where s G ?^S2- 

We give more details and additional examples for Tis^ in Appendix A. 3. 

3.2 Action of Dehn-twists and braid-moves 

The fusion basis states for a particular pants decomposition are simultaneous eigenstates of every 
Dehn-twist about a cylindrical segment in the decomposition. This is a main advantage of the 
fusion basis. 

For concreteness, consider a cylindrical segment with boundary labels i = (1, 1) and an anyonic 
fusion basis state \i,d) = |(1,1),t1) such that the segment is labeled = t and (i_ = 1 in the 
two labelings. Then the ribbon graph, before reducing to two dimensions, looks like 




(25) 



where we show only the cylindrical segment x [—1, 1], i.e., S2 x [—1,1], and not the portion of S 
to either side. By applying Eq. (16), a vacuum loop can be pulled out from a puncture within the 
shaded region, giving 




For clarity we have stopped drawing the outer boundary of the cylindrical segment x [—1,1]. 

The action of a Dehn-twist D on this cylindrical segment is to apply a full rotation. The Dehn- 
twist was defined in Section 2.3 as acting on S and ribbon graphs in E, but it is naturally extended 
to S X [—1, 1], and its action on ribbon graphs in E x [—1, 1] is easily seen to commute with the 
reduction to two dimensions S x {0}. The resulting state simplifies by moving the ribbon down to 
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S X {0}, pulling it over the vacuum loop with Eq. (16), and then removing the kink using Eq. (20): 




Thus D\{1,1),t1) = e~^'^*/^| (1, 1), rl). Had the two edge labels been both 1 or both r, then similar 
calculations would have shown \i, d) to be a +1 eigenstate of the Dehn-twist D. The eigenvalue 
is e"'"^'^*/^ when the edge labels are It. 

Now fix a pants segment of S, and let us compute the action of a braid move R in the anyonic 
fusion basis. Let d be a doubled anyonic fusion diagram, and say for example that the labels on 

TT TT 

the edges incident to the pants segment are given by '\</ . Then the corresponding ribbon graph 
in three dimensions is given by 

(28) 

where we have pulled a vacuum loop out from one of the boundaries of S. Applying R to this state 
yields 




R 




(29) 

Here we have again pulled edges across the vacuum line and applied Eq. (20) to remove the twists. 

Thus the action of R in this case is to swap the two portions of d beyond two of the boundaries 
of the pants segment, and add a phase. In fact, this would have been the action regardless of the 
doubled anyonic fusion diagram d, just with a different phase. The phase acquired corresponds to 
the i?-matrix of the doubled anyon model. Its general form will be derived in Section 6.5. 



3.3 Interpretation as a doubled Fibonacci anyon model 

The anyonic fusion basis can be interpreted in terms of the general theory of anyons, in particular 
of the doubled Fibonacci anyon model [FFN+OO]. This interpretation will justify our notation, but 
readers unfamiliar with anyonic theory may safely skip over it. 
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1x1 


= 1 


1 X T 


= T 


T X 1 


= T 


TXT 


= 1 + r 



The Fibonacci category Fib has two anyons, 1 and t. Each of the two labehngs of T defined 
in Section 3.1 is a Fibonacci anyon fusion diagram. An edge of T labeled 1 represents a vacuum 
anyon, while an edge labeled r represents a r anyon. A vacuum anyon can also be interpreted as 
the absence of a particle. The fusion-consistency constraint on a labeling requires that each internal 
vertex satisfy the fusion rules of Fib: 



(30) 



The final equation above should be interpreted to mean that fixing two of the incident labels at a 
vertex to be r and t, the third label can be either 1 or r. However, in the other three cases, fixing 
two of the incident labels at a vertex also fixes the third label. 

Notice that under both Dehn-twists and braid-moves the two ribbon graphs, in S x {1} and in 
S X { — 1}, from a doubled anyon fusion diagram transform independently of each other. Moreover, 
if these two ribbon graphs are swapped, then the phase of the eigenvalue of any Dehn-twist or 
braid-move is negated. 

Therefore, a doubled anyon fusion diagram d is the fusion diagram for the doubled Fibonacci 
category DFib = Fib ® Fib*. That is, DFib consists of two copies of the Fibonacci category with 
opposite chiralities. If we indicate the four possible pairs of labels for an edge by 11, Ir, rl or 
rr, then the allowed fusion rules are 



11x1t= It 
It X 11 = It 
Ir X It = 11 + It 



11 X TT = TT 

It X TT = TT + It 

t1 X TT = TT + t1 

TT X TT = 11 + It + t1 + TT 



11 X t1 
t1 X 11 
t1 X t1 



t1 
t1 
11 + t1 



(31) 



TT X 11 = TT 
TT X It = TT + t1 
TT X t1 = TT + It 



It X t1 
t1 X It 



TT 
TT 



More compactly, the above fusion rules can be written as 

(0102) X (6162) = (ai X 6i)(a2 x 62) 



(32) 



where aj,6j G {1,t}. They simply require that each coordinate separately be fusion-consistent. 

The fusion rules also allow us to determine the quantum dimensions of these particules, by 
regarding them as equations for the dimensions. The result is 



dii = 1 

dlT =drl=T 

d-TT = T + 1 



(33) 



where, recall, r = (1 + \/5)/2. 
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3.4 Changing the pants decomposition: The F-matrix 

To conclude this section, we would like to derive an expression for the basis change between anyonic 
fusion bases associated to different pants decompositions of the surface S. A pants decomposition 
of S corresponds to a rooted binary tree. A braid-twist swaps the subtrees at a pants segment, 
which allows for some basis changes. However, this only reorders subtrees, and does not change 
the parent-child structure of the tree. 

Any two pants decompositions having the same root puncture and the same left-to-right ordering 
of the leaf punctures can be related via a sequence of moves of the following type: 



(34) 



This move takes two pants segments connected by a cylindrical segment, and reconnects the middle 
subtree from one side of the lower pants segment to the other side. 

Let d be a doubled anyon fusion diagram for the pants decomposition on the left-hand side of 
Eq. (34). The change in basis to a superposition of doubled anyon fusion diagrams for the right- 
hand-side pants decomposition is most easily seen in the three-dimensional ribbon graph picture. 
From Eq. (3), we derive the two equivalent equations 



(35) 



By applying these rules separately to the ribbon graph in S x {1} and to the ribbon graph in 
S X {—1}, we obtain the anyon fusion basis change. As for Dehn-twists and braid-moves, the two 
ribbon graphs transform independently of each other. For example, we find that 




It It It 

It 
It It It 

1^ 



1 It It It i It It It 



1 It It It 



I It It It 



(36) 



since in the above cases the first labeling is always trivial, i.e., the ribbon graph in E x {1} is empty. 
To give an example of the basis change where both labelings are non-trivial, note that 



TT TT TT 




(37) 



In general anyonic theory, the matrix of this basis-change unitary is known as the F-matrix; in 
Fib Fib* , the F-matrix is the tensor product of the F- matrix for Fib with the F-matrix for Fib* . 
In fact, since the F-matrix for Fib has only real- valued entries, it equals the F-matrix for Fib*. 
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3.5 Recursive construction of the anyonic fusion basis states by gluing 

Let T^A = Sm+i and Sb = Sn+i be two spheres with n + 1 and m + 1 punctures, respectively. 
By gluing Ti^ to along a boundary component, one obtains the (m + n)-punctured sphere 
S = Tim+n- Pants decompositions for Tia and then combine to give a pants decomposition for 
S. Anyonic fusion basis states on S can thus be understood as the result of combining anyonic 
fusion basis states for T^a and S^- 

Indeed, for a labeling H. G {0, 1}™+"- of the boundary components of S, and a fusion-consistent 
doubled anyon labeling d of the pants decomposition, we have 

K,(i)s= «Ml4^)sJ^i^A>E^ ■ (38) 

fce{o,i} 

Here, l\ is the restriction of £ to the boundary components of Y^a^ additionally with label k assigned 
to the glued boundary; and (Ia is d restricted to the pants decomposition of T^a- The quantities 
£^ and dB are defined similarly. The coefficient a^^d depends on the anyon labels assigned to 
the connecting edge by d, according to 

1 if h+h- / TT 

i if h+h- = TT and A: = (39) 
A= if 6 I 6_ = TT and k = 1 

This can be derived by using the first identity of Eq. (35) to change the ribbon graph locally so 
that at most one ribbon crosses from Yja to S^. 

4 Realizing the Fibonacci ribbon-graph Hilbert space 

In this section, we define the Fibonacci surface code. Starting with a triangulation of a surface S, 
certain fixed boundary conditions and qubits placed on each edge of the triangulation, we describe 
commuting local projection operators such that their joint +1 eigenspace is the code space "Hi;. In 
Section 4.1, we define the local projection operators. These operators are essentially the same as 
those in the Levin- Wen Hamiltonian, except defined to take into account the boundary conditions. 
In subsequent subsections, we argue that the code space is and describe various properties of 
the code. 

4.1 The Fibonacci code 

Let S = S„ = S"^ \ (^4^ U • • • U A") be the n-punctured sphere, and let £ be a or 1 labeling of 
each hole A*. 

Let 7~ be a triangulation of S, and let T = (V, E) be the dual graph to T. 7" is a connected 
graph embedded in E. The components of S \ 7" are simply connected, and we refer to them as 
plaquettes. The vertices of T all have degree at most three. A vertex f in T can be univalent 
(respectively, bivalent) if the corresponding triangle in T has two edges (resp., one edge) along a 
boundary of S. It will be convenient to imagine boundary edges in the dual graph that correspond 
to the edges along the boundary in T. A boundary edge in the dual graph leaves a vertex v and 
crosses into some hole A^. If these boundary edges are included, then each vertex w in 7" is trivalent. 
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(a) (b) (c) 

Figure 2: The Levin- Wen model with a boundary, reahzing "H^: qubits sit on the sohd lattice edges. 
In (a) and (b), dashed edges represent "virtual" qubits fixed to 0. The thick, blue edge represents 
a virtual qubit fixed to the label i{p) of the marked boundary point p. In (c) is illustrated the 
gluing of anyon fusion basis states described in Section 4.2. The t virtual qubits are replaced by 
qubits, and prepared in states \k) \0)^^^^^\ for k G {0,1}, entangled with states on either side. 
Subsequently, all plaquette operators touching the boundary (shaded plaquettes) are applied. 

Now define a labeling on the boundary edges. For each boundary component, arbitrarily fix one 
boundary edge, and label that edge by the value of i on that component. Label all other boundary 
edges 0. We will refer to this new labeling also by i, cf. Figure 2. 

The Hilbert space Tij- = (C^)®'^' of our system is obtained by placing a qubit on each edge e. 
Note that boundary edges are not given qubits. We use orthonormal bases {|0)e, |l)e} for the local 
Hilbert spaces. As in Section 2.2, we will use the convention that a 1 indicates the presence of a 
ribbon along the edge, while a indicates the absence of a ribbon. 

The Fibonacci code is a subspace of (C^)®^, defined as the simultaneous +1 eigenspace of a set 
of commuting projections. There are two different kinds of projections: 

• For every vertex v of the graph, there is a vertex projection Qy. Qy depends only on the 
three edges incident to v, and on those three edges is the diagonal operator 

Qv= Yl . (40) 

i,j,k£{0,l}: 

That is, Qv imposes that the vertex u is a legal ribbon graph vertex without dead-ends. If v 
is connected to one or two boundary edges, then is defined as in Eq. (40), except with the 
boundary edges fixed by the labeling L For example, if there is one boundary edge e, then 

is given by Ei,j:i+j+£(e)^i NiX^il on the other two edges. 
Since all of the vertex projections are diagonal in the same basis, they commute with one 
another. Let P be the number of plaquettes in T and let A = S„+p be the surface obtained by 
placing a puncture in the center of each plaquette in S. Then the simultaneous +1 eigenspace 
of the vertex projections is T-L^^'^ ^ , the space of valid ribbon graphs on A with open boundary 
conditions on the new punctures. The computational basis for the qubits corresponds to the 
computational basis for defined in Section 2.2. We call this space the physical subspace 

of n^. 



16 



In more detail, we regard every physical configuration as a ribbon graph on A, using the 
embedding of T into A. Because of the punctures V, no two different computational basis 
states or physical configurations are identified with the same ribbon graph. This shows that 
physical subspace is indeed isomorphic to H.)^ . 

• For every plaquette p of the graph, there is a plaquette projection Bp. Bp is supported 

on 1-6^^ ^ and acts on it by adding a vacuum loop, divided by -v/T+T^, around the puncture 
in the middle of p. Its action on (C^)®'^ can be obtained by reducing the resulting ribbon 
graph back to the computational basis. This definition of the plaquette projection Bp in 
terms of a vacuum loop is given in the appendix of [LW05b] (see Section 6.6 for more details). 

Bp depends on all the edges and boundary edges incident to the vertices around the pla- 
quette p. However, it only possibly changes the qubits on the edges circling the plaquette, 
i.e., it is a controlled projection, controlled by the edges with exactly one endpoint on the 
plaquette boundary. 

Since the vacuum loop added by Bp is separated from other vertices or plaquettes, the pla- 
quette operators commute with each other and with the vertex operators. 

By Eq. (16), the plaquette operator Bp has the effect of removing the puncture from plaquette p, 
allowing ribbons to be pulled across it. Therefore, the simultaneous +1 eigenspace of all of the 
vertex and plaquette operators is isomorphic to T-L^ (see Lemma 4.1 below). Thus we have defined a 
quantum error-correcting code (QECC) based on ribbon graphs. It encodes Ti.^ inside of (C^)®''^'. 
We denote the code space by c 71;^. Here the superscript gs stands for ground space., because 
T-d^^ = is the ground space for the Levin- Wen Hamiltonian 

Hi = - Yl Bp- Y . (41) 

plaquettes p vertices v 

To describe the isomorphism between T-L^ and the simultaneous +1 eigenspace 7^^'*^ ^ of all 
plaquette- and vertex-operators on T-Lj- = (C^)'^l'^l in more detail, let B be the projection 

B = llBp. (42) 

p 

Then the following lemma, proved in Appendix B, describes how ribbon graphs are mapped to the 
qubits on T. 

Lemma 4.1. Given a state |^'e) G 'Hy.j deform the ribbon graph (arbitrarily) to avoid all the points 
in V. The state can then be regarded as an element £ Ti^^'^ ^ . Then -B|<I>a) does not depend 
on the initial deformation (i.e., \^a))- The map 



K:W^^ (43) 



f 

|^e)^B|<5a) , (44) 
is an isomorphism (preserving the inner product). 
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Lemma 4.1 allows us essentially to forget about the graph T, and work completely in the off- 
lattice picture Ti^. For any state in H^, the corresponding codeword, or ground state, |^'^) € Ti^jF^ 

is obtained by mapping the off-lattice ribbon graph to the graph T, getting a state I'&a))) and then 
applying the projection B. 

It will often be convenient to use the state |$a) to represent the ground state l^*^) = ^I'&a)- 
Of course the projection B is not one-to-one. However, one can argue that B\^^) = B\^'^) if and 
only if |$a) and |$^) are equivalent under certain discrete lattice versions of the local relations in 
Eqs. (l)-(3) (see [LW05b]). 



4.2 Anyonic fusion basis states and gluing 

Here we describe the discrete analog of the gluing property explained in Section 3.5. Let T be a 
triangulation of a surface S that is obtained by gluing to along a closed curve 7. Assume that 
7 is part of the triangulation T. Let Ta and Tb be the restricted triangulations of T to and T,b, 
respectively, and Ta and Tb their dual graphs. Figure 2 shows these dual graphs schematically; 
edges intersected by 7 are shown as "virtual edges" attached to the boundaries of and E^, 
respectively. Pick an edge e from among the t virtual edges. 

The images \£, d)j- = A|£, d) G of the anyonic fusion basis vectors \l, d) G satisfy 

essentially the same properties as their off-lattice counterparts. Observe that plaquette operators 
within Tja and do not affect qubits on the virtual edges, and can be applied before the plaquette 
operators ^(7) = np-pn7^0 -^p "-"^ plaquettes intersecting 7, by commutativity. If we apply the 
isomorphism A to a state decomposed as in (38), the result can therefore be written as 

\i,d)^ = B{^) aMKidA>rj4,rfB>tJfc)e|0>^(*-^) , (45) 

fee{o,i} 

illustrated in Figure 2(c). Here, the states 1^^,^^)^, |^B'^B)fB ^'^^ anyonic fusion basis 
states associated with Hamiltonians -ff^, HJ^, respectively. Thus \l,d);^ results from applying 

projections along 7 to a certain superposition of product states on ® ^Ta i'^'^)^^- 

Consider now the case t = 1. Then -6(7) = Bq for a single plaquette q. The following lemma 
states that the qubit on edge e directly reveals some information about the anyon labels 6+6_ 
assigned to the connecting edge between T,a and by the doubled fusion diagram d. 

Lemma 4.2. Let T , Ta, Tb be graphs as described above with t = 1, and let e be the connecting 
edge. Then the anyonic fusion basis state \i, d)j- is given by 



r 



\(.%dA)f^\fB,dB)f^\^)e i/6+6_ = ll 

K:dA)f^\l'B,dB)f^\l)e i/6+6_ e{lT,rl} (46) 

>l^^,rfA)tJ^^,dB)tJO)e + ^|^il,dA)tj4,tiij)tJl)e ^/ M 

A proof of this statement is given in Appendix B. 



TT 



4.3 Relating ground spaces for different triangulations 

In preparation for presenting a scheme for computing on codewords, in Section 5, let us conclude 
this section by deriving how the code space changes under local changes to the triangulation T. 
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Figure 3: Change of triangulation T by flipping edge e (left), and the effect on the dual lattice T. 



Let e be an edge in T that does not go along the boundary of S. Let T' be a triangulation 
obtained from T by reconnecting the edge e as in Figure 3. In the dual graph T, this corresponds 
to an operation on up to five edges to construct T' ■ 

Since e does not go along the boundary of S, we can use the same boundary conditions £ for 
T' as for T. Moreover, the dual graph T' has the same number of edges \E\ and punctures P, and 
can be used to define a punctured surface A'. Thus as before we can define Hilbert spaces T-L^jt,, 

y.^, a Hamiltonian Hi^^, and the associated subspaces H^^f ^ and T-L^^\ 

To study the relationship between the two ground spaces n'ii^' and ^^f , we define a linear 




(47) 



in the computational bases. Here, a solid line represents the state whereas a dotted line 
represents the state |0). We can easily extend this to a unitary Fg : Tij- — )• (To make sense 

of unitarity, use the obvious isomorphism between the two spaces, each isomorphic to (C'^)®'^'.) 
We will henceforth refer to this as an F-move; depending on which is more convenient, we will use 
either e or the dual edge e to label this move. 

It can be shown [KRV09] that maps the ground space 'H^^'* isomorphically to Ji^^^ More- 
over, this isomorphism is compatible with the representation of ground states from Lemma 4.1 in 
the following sense: 
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Lemma 4.3 ([KRV09]). The following diagram commutes, where B' is the product of the plaquette 
terms for T' , as in Eq. (4^): 



iL A iLa 



T 



(48) 



^(£,0^) B' _ ^^gs 
rt^l tt;: 



T' 



5 Computation with doubled Fibonacci anyons in Ti^, 

In this section, we explain how to operate on the code: how to encode qubits, prepare and measure 
states, and apply unitary operators on codewords. 

Universality results for Fibonacci anyons were obtained by Freedman et al. [FLW02, LW05a], 
using an encoding of a logical qubit into three-anyon states: 

T T T T T T 



This encoding was also used by Bonesteel et al. [BHZS05], who showed how to efficiently approx- 
imate any gate in terms of a sequence of braids, using the Solovay-Kitaev construction [KSV99]. 
This encoding, while efficient, has the physical disadvantage that the total anyonic charge of the 
code-state is non-trivial. Hormozi et al. [HZBS07] used an alternative encoding, with four r anyons 
for each qubit: 



T T T T 

T T X[ T T 




|0)^ (50) 

1 1 

As explained in [FLW02, BHZS05, HZBS07] (cf. [Kit03, Pre04, NSS+08]), an arbitrary n-qubit 
circuit can be efficiently approximated if we can 

• Create the encoded state |0)®". In both encodings (49) and (50), this amounts to preparing 
the state of n pairs of r anyons, each pair fusing to the trivial particle 1. 

• Execute arbitrary braids. 

• Measure the anyon label {1, r} of a fixed edge in a fixed fusion tree basis. 

We will use the "doubled" versions 11 and rl instead of 1 and r. The other anyons in the 
doubled theory. It and rr, will not be used for computation. 



5.1 State preparation and measurement 

A method for preparing the (unique) ground state of the Levin- Wen Hamiltonian on a sphere has 
been given in [KRV09] . The essential idea is to apply a version of discretized surgery, which easily 
generalizes. 

Assume that we have a codeword l^*) for the Fibonacci code on a surface S, and a codeword 
1$) on the surface T. Then [KRV09] gives a way of cutting a hole inside a plaquette of S, and a 



20 



Figure 4: The procedure cut{'y) isolates a region enclosed by a closed curve 7, by making 
F- moves in sequence along the edges ei, 62, . . . , e„, colored red, counterclockwise along 7 on the 
dual graph T. The result is a tadpole-like structure as in Lemma 4.2, with a single edge e and 
plaquette q between the two regions. The inverse operation is glue{'y). 

hole inside a plaquette of T, and gluing these holes together. Here we discuss the details needed 
to create specific codewords on the n-punctured sphere. 

As explained in Section 4.2, the gluing of codewords associated with two surfaces T,a and T,b 
takes a particularly simple form for graphs with tadpole-like structures. To make use of this fact, 
we change the triangulation to obtain such structures: Figure 4 shows how a region T,a bounded 
by a closed curve 7 can be isolated from the rest by a procedure cut('y) consisting of F-moves. The 
inverse procedure glue{'j) integrates a tadpole-like structure into a regular lattice. 

Procedure cut{j): Let ei, 62, . . . , G T be the sequence of edges in counterclockwise order con- 
stituting a closed path on T along 7. For i from 1 up to n, apply Fe-. 

Procedure glue{'y): Let e[,e'2, ■ ■ ■ ,e'^ be the edges in the deformed lattice corresponding to the 
sequence ei, 62, • • • , in the original lattice. For i from n down to 1, apply F„' . 

i 

Figure 5 shows the intermediate steps in an example of applying cut{'y). Note that the F-moves 

constituting cut{'y) do not commute with each other. Note also that these F- moves lead to a 

degenerate triangulation. By adding more F-moves, this procedure can be modified to keep the 

degree of the bounding plaquette small. 

Our goal is to create 2n anyons of type rl, which fuse to 11 in pairs. This can be achieved by 

rl rl 

n times cutting out a circular disk and replacing it by a three-punctured sphere with state N<;|^ . 

11 

To describe this in more detail, let 7" be a triangulation of S, a surface obtained by gluing Tia 
to S^. Let d-A and be fusion diagrams with each root edge carrying the trivial label 11, and 
consider the state {lAis, dAdB)f- 

The necessary steps for cutting out T,a and replacing it by a different surface Y^a' ai'e especially 
simple when T has the tadpole-like form considered in Lemma 4.2; we merely have to trace out the 
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i = 4 



i = 5 



i = 21 



Figure 5: The successive figures show the intermediate steps of the cut{'y) procedure, after applying 
Ff,- . The thick, red edg es are ei, . . . , e^' The overall effect of the procedure is shown in Figure 4. 



qubits on Ta, and replace them by a state |^a'; f^A')s^/ • Indeed, by Lemma 4.2, we have 

\iAiB,dAdB)f = \e%dA)fJl%,dB)fJO)e 

S|^°S,dB)tJO)e 

-f^''\l%,dA')fyB,dB)fJO)e 
= \iA'^B,dA'dB)f, ■ 



(51) 



To give T the tadpole form, choose any closed curve 7 not enclosing any holes in the surface — this 
ensures that the corresponding fusion diagrams dA and dB have trivial total charge. Then use the 
cut{'y) and glue{'y) procedures as in the following circuit: 

Tb 



B 



dAdB) 





cut{'y) 


e 



















(52) 



\iA',dA')f^^ 

'A' 

In our application, the triangulation Ta' can be chosen to be a coarse triangulation of S3, 
such that the state \£A'TdA').f- corresponding to a rl anyon pairs can be created by a projective 

'A' 

measurement on a small number of qubits. 
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According to Lemma 4.2, the qubit on the tail of a tadpole is or 1 according to whether the 
anyon label of the corresponding pant segment is 11 or rl, respectively. This leads to the following 
circuit for measuring the anyon across a particular pant segment: 

Tb 





cut{'y) 


e 














Ta 



(53) 



Implementing these procedures in the presence of noise, with periodic error correction, both 
measurements and preparations expose the encoded information. Beneath a threshold, a constant 
measurement error rate can be reduced by taking the majority of multiple measurement outcomes. 
Errors in preparation can be efficiently dealt with using the composite Fibonacci anyon distillation 
scheme of [K6n09], in order to concentrate the entropy into specific regions using reversible gates. 



5.2 Implementing diffeomorphisms on S by changes of triangulation 

Recall that the mapping class group of S acts on by deforming embedded ribbon graphs. If 7" is 
a triangulation of S, such that the dual graph T has N edges, then Hj- = (C^)®^. In this section, 

we show how to implement the mapping class group's action on 7^^^* = — both Dehn-twists 

and braid-moves — by local unitary operators on (C^)*^^. As before, the implementation is based 
on changes of triangulation using the operators F(, from Eq. (47). 

Consider first a Dehn-twist, specified by a simple closed curve 7, as in Eq. (12). Assume that 
7 is supported on edges of the triangulation T, and let T^l and S^j be the triangulated surfaces on 
either side of 7. The Dehn-twist operator -D(7) is obtained by applying F-moves along edges of 
one of these surfaces, more precisely by the following procedure. 

Dehn-twist D{-j): Let 7 consist of n edges. Repeat the following for rounds r = 1, . . . , 2{n — 1): 

In round r, consider the triangulation T"*""^ constructed after r — 1 rounds. = T is the 
initial triangulation. Let Tq, . . . , T„_i be the triangles bordering 7 in S/j, in counterclockwise 
order. Assume that each triangle has exactly two vertices on 7; if true in the first round, then 
this will hold always. Let Vi be the vertex shared between Tj and Tj+i (mod n) ) let Wi be the 
vertex of Tj not in 7, and let be the edge {vi, Wi). For i from 1 to n, apply the maps Fg. to 
obtain the triangulation 1"^ . 

Note that the F-moves within each round commute, so can be applied in an arbitrary order. 
Figure 6 shows an example of this procedure. A clockwise Dehn-twist can be defined similarly. 

With the isomorphism A described in Lemma 4.1, and the explicit form (47) of the maps 
and Lemma 4.3, it is straightforward to see that the operator .0(7) implements a ribbon-graph 
Dehn-twist on T-Ly. — T~i^, as studied in Section 2.3. Figures 7 and 8 show an example, tracking 

through F- moves states |<I>) G T-i^^^^ that represent ground states. 

A half-twist twists along a curve 7 by vr instead of 2tt. Define an operator ^ D{'^) on (C^)®^ to 
execute only rounds r = 1, . . . , n — 1 of D('y). The braid move of Eq. (11) can then be implemented 
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(a) r° (b) (c) 

Figure 6: From left to right, the triangulations T^, and obtained while implementing a 
Dehn-twist about the inner curve 7. The portion of the surface within 7 is not shown. The 
thick, red edges are the locations of i^'-moves in the next round, i.e., the edges {ej}. 

by three half Dehn- twists about the curves 7a, 7b and 7c given by 




Indeed, we have 




(55) 



Therefore both braids and Dehn-twists on the ribbon-graph Hilbert space can be implemented 
using F- moves. 

6 Quantum computation with Turaev-Viro codes from general 
anyon models 

Although for concreteness and to minimize notational complexity we have focused on the doubled 
Fibonacci anyon model, the calculations and methods generalize. In this section and Appendix B, 
we will go over the construction of surface codes for other doubled anyon models. Section 6.3 
begins by collecting the data needed to specify an anyon model. Section 6.4 introduces the ribbon 
graph Hilbert space associated with a surface S and includes the derivations of a number of 
useful local identities. Section 6.5 introduces the doubled fusion basis states by considering ribbon 
graphs in S x [—1,1]. Section 6.6 describes how the ribbon graph Hilbert space is realized as 
the ground space of the Levin- Wen Hamiltonian. 
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Figure 7: This figure shows for various i how i F- moves in the first round (r = 1) of the Dehn- twist 
D(7) affect a ribbon graph |$) G ^^'^ ^ embedded in the dual graph T. 




Figure 8: This figure continues Figure 7 for larger i, in order to illustrate that D^-j) indeed cor- 
responds to the deformation of an off-lattice ribbon graph by a Dehn-twist. The half-Dehn-twist 
Y^D(7) is implemented by i = 56 F-moves. After the 107-th F-move, we switch from |$) G Ti^^'^ ^ 
to a different representative \^') G Ti^^'^ ^ for convenience (note that = i?|^>')). 
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6.1 Background on categories 



While we only rely on basic definitions of "anyon models" for the generalization discussed below, 
some category-theoretic terminology is required to set our work in context. A thorough discussion of 
these concepts is beyond the scope of this introductory article. We restrict ourselves to a summary 
and refer the interested reader to the literature for detailed definitions. 

Roughly, a spherical category as introduced by Barrett and Westbury [BW99] is a set of data 
which gives rise to an isotopy invariant of trivalent directed labeled planar graphs embedded in the 
two-sphere. A corresponding "coherence theorem" [BW99] guarantees consistency (of local rules) 
under a finite set of conditions on the data. More specific properties are needed in our context: 
a unitary braided fusion category (see, e.g., the appendix of [Kit06]) additionally has a notion of 
braiding, and is equivalently referred to as a unitary ribbon category. The reason for this terminology 
is that the braiding notion allows to make sense of crossing edges (i.e., non-planar graphs) and also 
allows to introduce a notion of twisting. This leads to an invariant of ribbon tangles, as studied by 
Reshetikhin and Turaev [RT90] . The coherence theorem in this context is more commonly referred 
to as Mac Lane's theorem [Mac98]. Finally, a unitary modular tensor category is a unitary ribbon 
category with a non-singular S'-matrix, as defined below; the terminology here stems from the fact 
that this give rise to a (projective) representation of the modular group. 

There are multiple ways of obtaining new categories from old ones. The categorical double or 
Drinfeld centre DC (see [MiigOSb] or [Kas95, Section XIII. 4]) of a (not necessarily braided) fusion 
category C is always a braided category. Moreover, if C is spherical, then DC is modular, as shown 
by Miiger [MiigOSb]. The category DC has a particularly simple structure if C is itself already 
modular: in this case, DC = C C* is isomorphic to the direct product of C and the conjugate 
category C* obtained by complex conjugation of all data (see below). We formulate all our results 
for this case. 

6.2 Quantum codes from categories 

A Turaev- Viro TQFT can be defined using any spherical category C. In particular, the arguments 
given below in Section 7.2 imply the following: For any finite unitary spherical category C, and 
for any triangulated surface S, with or without boundary, there is a corresponding quantum error- 
correcting code. This code has qudits assigned to the edges, where d is equal to the number of 
simple objects in C. Unless C is multiplicity-free, the triangles also need qudits. 

For certain categories C, this code can be expressed as the ground state of a Levin- Wen local 
stabilizer Hamiltonian [LW05b]. Here we focus on the case where C is modular, and show that the 
anyons of the code are from the category C C*. A more general theorem [TurlO], proved in the 
modular case in [Wal91, Tur94], asserts that if C is any finite semisimple spherical category, then 
the Turaev- Viro or Barrett- Westbury TQFT is the same as the Reshetikhin- Turaev TQFT for the 
quantum double DC (see e.g., [MiigOSa]). 

6.3 Parameters of an anyon model 

Consider an anyon model described by a tensor category C. Such a model has particle types 
{l,i,i* , j, j* , . . .}, where * denotes charge conjugation and 1 = 1* is the trivial particle. Let 



di = di* be the quantum dimension of particle i (note that di = 1), and let V = yj^i^i be the 
total quantum dimension. Let 5abc* = 1 if the fusion space V^f^ is at least one-dimensional, and 
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otherwise. Assume for simplicity that the anyon model has no fusion multiplicities. (Our results 
directly generalize to anyon models with fusion multiplicities.) These fusion rules satisfy 

y^,Sijm*Smkl* =y^,Sjkm*Siml* > (56) 
m m 

which expresses associativity of fusion. The quantum dimensions satisfy 

didj = ^ Sijkdk . (57) 
k 

It is further necessary to specify the operations for fusing and braiding anyons. The fusion 
tensor, or i^'-move, for changing fusion bases has entries 

i j k ^ , . . , i j k 

It satisfies the five properties, for all s, 



physicality: 




^ki^^ijm^Mm* 




i^jkn* 


pentagon identity: 




rpmlq rpjip* rpjsn 
^kpn -^™"«-^£fcr 


— pjip* pi"*iq* 
q*kr mis 


unitarity: 








_ rpi*j*m* 
— -f 




tetrahedral symmetry: 




Tpijm 
^kin 


_ rpjim 
~ ^ikn* 


Tplkm* _ 

jin 


_ zpimj 
" ^k'ni 


normalization: 






Tpii*l 




The braiding operations, called R- 


-moves, are f 


;iven by 








J ' 

Y 






i 

ii^k'*r 


i ■ 


k 


k 




k 




k 
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(59) 



(61) 



(60) 

k 

The entries satisfy \R^^ \ = 1 and the hexagon identities: 

jjki T?k* i* m Tjkj \ ^ Tpi*k''m jjkn Tpj*i*n 

n 

Tfik TTikim* -pjk \ ^ rpikm* jjnk injin* 

^m^e*jg*^g — 2^ ^ i* jn* -"-t ^ E*kg* 
n 

Combined with F-moves, i?-moves allow for the resolution of crossings, via, e.g., 

Xr^{i'"-''^'X,- '''' 

They also allow for defining the topological phases and the topological ^-matrix. The topological 
phases are given by 9i = {R\ *)* and allow a ribbon to be untwisted: 



Ci=ei'-i. (63) 
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The topological 5-matrix is defined by 



= ^ ^QD^ ■ (64) 

This matrix is symmetric and unitary; it satisfies 

Sij = Sji = Si*j* = Sj*i*, Sil = di/T> . (65) 

For an anyon model with particles {1+, a+, 6+, . . .} described by a tensor category C+, we can 
define a dual category C_ with anyons a_, . . .}. Here a_ is the same as a+, but with opposite 
chirality, i.e., the topological phase is given by 9a_ = 0*^. Furthermore, the i?-matrix is replaced 

by its adjoint, which for the multiplicity-free case considered here boils down to RaZ'^~ = {Ra'^'^^)* ■ 
The double of C = C+ is the category DC = C+(8>C_; its particles are {a+(8>6_ | a+ G C+,6_ G C_}, 
and the braiding and fusion matrices are tensor products of those of C+ and C_ . (Note: Here we 
used the assumption that C is modular. For a general category C, the double DC need not be 
isomorphic to C+ C_.) 



6.4 The ribbon graph Hilbert space "H^ 

Using the fusion rules 6ijk, the quantum dimensions di and the F-tensor Fl^^ , we can define the 
ribbon graph Hilbert space 1-Ly. associated with a surface S as in Section 2.1. A (colored) ribbon 
graph is a graph with labeled directed edges embedded into S, with vertices of degree two and 
three in the interior, and degree-one vertices at the marked boundary points of S. Switching the 
direction of an edge corresponds to conjugating the anyon label, and edges assigned the trivial 
label 1 can be removed or added from the picture according to 




Valid ribbon graphs are those that satisfy the fusion rules at every vertex, i.e., having 5ijk = 1 
for every vertex with incoming edges labeled i, j and k. (For a vertex of degree two, this rule 
is adapted by adding an edge with label 1.) The Hilbert space T-Ly, is the space of formal linear 
combinations of such ribbon graphs, modulo the local relations 




Compare to Eqs. (l)-(3). 

The Hilbert space decomposes into spaces T-L^ indexed by labelings £ of the marked boundary 
points on the boundary components of S, as in Eq. (7). is the subspace of ribbon graphs with. 
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for each boundary point p, an edge carrying label £{p) leaving p. An orthonormal basis of can 
be constructed using ribbon graphs living on the dual graph of a triangulation of S, as explained 



in Section 2.1. For example, two orthonormal bases of the space 'Hj,^ 
by the ribbon graphs 



{k*,e) 



on the annulus are given 



B' 



Ski*j* / 0, 6ij£* / 



(70) 



The inner product can be defined using any such basis; the result is independent of the chosen basis 
because different triangulations of S with the same number of edges are related by edge flips and 
i^-moves are unitary. 

6.5 Ribbon graphs on S x [— 1,1] and fusion basis states of the doubled model 

The anyonic fusion basis states for the doubled theory are defined, as in Section 3.1, by reducing 
"doubled" ribbon graphs in the manifold S x [—1, 1] down to the surface S. 

In the following, we use undirected, unlabeled, dashed lines to denote the superposition of all 
anyon types, weighting i by di/V: 



(71) 



These vacuum lines generalize Eq. (15). Lemma A.l below summarizes their main properties. 

Consider a fusion diagram describing a state of n anyons of a doubled theory DC = C+ . 
This is a tree G with n + 1 leaves at the punctures of the (n + l)-punctured sphere S = Sn+i, 
and with edge labels of the form ® j_. Consider the manifold M = S x [—1, 1], and take two 
copies of G, G± = G x {±1} C S x {±1}. For an edge e G G labeled i+ (8) j- in the original 
fusion diagram, label the corresponding edges e+ and e_ in G+ and G_ by i and j, respectively. 
In addition to "doubling" the edges in this fashion, add vacuum lines, from Eq. (71), embedded in 
S X {0} around n of the punctures. 

For concreteness, consider the case of a single anyon i+ (8) j_, corresponding to the ribbon graph 




(72) 



on T,2 X [—1, 1]. In this figure, the shaded area is the inner component of {dT,2) x [—1, 1]. 

To see that this is consistent with an (8 j_ anyon, observe first that it has the required 
idempotency property, since stacking gives the same diagram up to a constant factor: 



V 



(73) 
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according to the "doubling" property (91) of vacuum lines stated in Lemma A.l. We can also 
verify that it carries the correct topological phase 9i_^6j_ = di^j, by computing the action of a 
Dehn- twist D('j): 




Here, as in Eq. (27), we have used Eq. (90) to pull the ribbons over the vacuum line, and then 
removed the kinks using the definition of the topological phase, Eq. (63). 

Next, let us check that this realization reproduces the correct i?-matrix, by considering a braid: 




Using Eq. (60), we find that, indeed. 

It is easy to see that this definition of "doubled" ribbon graph on S x [—1, 1] also gives rise to 
the correct F- matrix of the doubled theory DC, since the tensor factors behave independently. 

Note that this realization of the doubled theory by ribbon graphs on S x [—1,1] is only an 
intermediate step; we are ultimately interested in ribbon graphs on S. Ribbon graphs on S x [—1, 1] 
can be reduced to S. First, pick a label i{p) for every boundary point p G dT, fusion-consistent 
with the anyon labels i+ at {p, 1) and i_ at (p, —1), respectively. Then, connect up these two ends, 
and attach an edge with label £{p) ending at (p, 0). Clearly, the resulting ribbon graphs have only 
one edge ending at every hole; in particular, we can project the string-net onto S, and then resolve 
crossings by use of identity (62). This eventually gives a ribbon graph on S, an anyonic fusion basis 
state. 

For example, for S2 we obtain from Eq. (72) the ribbon graphs 




provided dki*j* / and 6ij£* / 0. Such ribbon graphs are explicitly reduced to elements of 'He2 
for the Fibonacci theory in Appendix A. 3. By Lemma A. 2 below, the ribbon graphs /c^'^(«, j) are 
idempotents under stacking. To give another example, the described recipe gives basis vectors for 



30 



H.Y.i of the form 




(78) 



Appendix A. 2 proves that the anyonic fusion basis states form a complete, orthonormal basis. 



6.6 The Turaev-Viro code as the ground space of the Levin- Wen Hamiltonian 

In this section, we present Levin and Wen's lattice Hamiltonians [LW05b] realizing the space T-L^ of 
(colored) ribbon graphs for a general anyon model with no fusion multiplicities, as in Section 6.3. 
This generalizes the construction of Section 4 for the Fibonacci model. 

As in Section 4, let T be a triangulation of S, and let T = iy,E) be the dual graph of T. 
Let rtf = be the Hilbert space obtained by associating a A^-dimensional qudit with 

orthonormal basis corresponding to the different anyon labels to every edge. Computational 
basis vectors will be represented by directing the edges of T and labeling them by anyon labels; 
reversing the direction of an edge corresponds to conjugating the label. Labelings I of boundary 
points p G 9S by anyon labels will be extended to virtual boundary edges attached to T as discussed 
in Section 4; that is, we associate the trivial label 1 to all but one edge that carries label ^{p). 

On T-L;f, we define vertex and plaquette projections as follows: 

• For every vertex v of T, there is a vertex projection Q^. This operator depends only on the 
edges incident on v, and enforces the fusion rules a.i v. It is given by 

Qv = ^^5ijk\ijk){ijk\ . (79) 

As before, these are commuting projections. Their simultaneous +l-eigenspace is (isomor- 
phic) to the space H}^ of ribbon graphs on a surface A. Here A = S„+p is obtained by 
placing a puncture into each plaquette of T embedded in S^. Let 1^ be the set of these new 
punctures. The ribbon graphs in have open boundary conditions on punctures p £ V, 

and are in one-to-one correspodence to the computational basis for the qudits. 

• For each plaquette p of T, there is a plaquette projection Bp. Bp is supported on the sub- 

space T-L\' , and acts by adding a vacuum loop, divided by total quantum dimension V. By 
reducing the resulting ribbon graph, its action on computational basis vectors is given by 

Bp = ^^diOl where 

i 

Here we have assumed that p has r boundary edges, and have identified jo = > and = k^. 
The plaquette operators commute with each other and with the vertex operators. 




fci h 



•rrivjvlv-i 
k 1 k 




(80) 
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The Levin-Wen Hamiltonian on on T-L^ is the sum of these operators, as in Eq. (41). 
According to Lemma 4.1 (which holds for the more general definitions given in this section), its 
ground space c T-i^^, i.e., the simultaneous +l-eigenspace of all plaquette and vertex op- 

erators, is isomorphic to the ribbon graph space T-Ly,- Quantum computations can therefore be 
performed in a similar manner as before using F-moves, i.e., five-qudit gates. A generalized version 
of Lemma 4.2 allowing for state preparation by gluing and measurement of topological charge is 
given as Lemma B.2 in the appendix. 



7 The Turaev-Viro invariant and the Turaev-Viro code 

In Section 6.6, we have discussed how the ribbon graph Hilbert space T-Ly, can be realized as the 
ground space of the Levin- Wen qudit lattice Hamiltonian. Here we give a more direct formulation 
that justifies calling this ground space the Turaev-Viro code. This approach historically preceeds 
the Hamiltonian construction and explains its origin. It is the extension of the Turaev-Viro invariant 
of 3-manifold to a TQFT [TV92]. 

We proceed as follows. In Section 7.1, we define the Turaev-Viro invariant [TV92] in the 
general category-theoretic formulation due to Barrett and Westbury [BW96]. Section 7.2 defines 
the Turaev-Viro code for a surface E, based on the Turaev-Viro invariant for the thickened surface 
S X [—1,1]. Section 7.3 shows that this definition corresponds with the Levin- Wen Hamiltonian's 
ground space, by deriving the Levin- Wen Hamiltonian; applying a plaquette operator corresponds 
to attaching a "blister" to the surface. 

Appendix C gives an alternative definition of the Turaev-Viro invariant, due to Kontsevich, 
Crane and Kohno [Kon88, Cra91, Koh92], and sketches how the anyonic fusion basis state decom- 
position of T-Ly. can be used to argue that the different definitions are equivalent. This background is 
useful for the proof in [AJKRIO] that approximating the Turaev-Viro invariant is a BQP-complete 
problem. 



7.1 Definition of the Turaev-Viro invariant 

In this section, we present the definition of the Turaev-Viro invariant TVc(M) of an oriented 3- 
manifold M. For simplicity, we assume that the F tensor of the category C satisfies tetrahedral 
symmetry, though this is not necessary (see [BW96]). The manifold M may also have a boundary, 
in which case the Turaev-Viro invariant TVc(M, x) depends on X; a labeling by anyon types of the 
oriented edges of a triangulation of the boundary. 

Triangulate M compatibly with the possible boundary triangulation. For each tetrahedron 
place a total order <t on its four vertices, such that the orders are consistent across tetrahedra, 
i.e., if t and t' share an edge {v^w}, v <t w when v <t' w. For example, such orderings may be 
obtained from a single ordering of all the vertices in the triangulation. 

For a labeling (j) of the edges of the triangulation by anyon types from C, and for a tetrahedron t, 
define a scalar gf as follows. First, we consider the "standard" tetrahedron, with vertices ordered 
vq < ■■• < vs and labels i,j,k,£,m,n of the six edges as shown in Figure 9(a). Here we have 
oriented the edges according to the vertex orders. Define the tensor 



I J m 
kin 



pi jm 

. (81) 
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V3 
(a) 



^2 
(b) 



(c) 



Figure 9: A standard tetrahedron labeling and edge orientations (a), and two other labelings. 

Now for a tetrahedron t, the value gf is obtained from this tensor by arbitrarily aligning t to the 
standard tetrahedron, and then conjugating a label if the orientation of the edge does not agree 
with the standard orientation. For example, if t is the tetrahedron of Figure 9(b), then 



9t 



I J m 
k* t n* 



(82) 



The tetrahedral symmetry of F implies that gf is well-defined. 
Then set 



TVc(M,x)=p-^i^-w^-^i y: n n v^0w n at, (83) 

labelings <j>- edges edges tetraliedra t 

<l}\aM=X eeM\dM eedM 

where \Vm\ is the number of vertices in the triangulation, and iVgjv/l is the number of vertices on 
the boundary. The sum is over all labelings (j) that agree with x on the boundary after conjugating 
labels so edge orientations on the boundary triangulation agree with the orientations from the full 
triangulation. It can be shown that the quantity TVc(M, x) does not depend on the choice of 
triangulation in the interior of M. 



7.2 The Turaev-Viro code 

Let S be a triangulated two-dimensional surface. For simplicity, assume that T, lacks a boundary. 
Consider the 3-manifold M = T, x [—1, 1]. Triangulate its boundary dM = S x {—1,1} with two 
copies of the triangulation of S. Then a labeling x of this triangulation consists of a pair (x+, X-) 
of labelings of the triangulation of S. In particular, we can interpret the quantity TVc(M, x) as 
the matrix element of an operator 

Wc(Sx [-1,1])= Yl TVc(M,x)|x+>(X-l • (84) 

x=(x+.x-) 

This operator acts on the Hilbert space (C™)®l-^l of labelings of the \E\ edges of the triangulation 
of S by the m anyon labels of C. It only depends on the boundary triangulation, and not on a 
triangulation of M. 

Define the Turaev-Viro code on S as the range of TVc(S x [—1, 1]). Since gluing together two 
copies of S X [—1, 1] along a boundary component gives again T, x [—1, 1], and since the gluing of 
manifolds corresponds to the concatenation of the linear maps from Eq. (84), TVc(S x [—1,1]) is 
a projection. 

We will next show that the Turaev-Viro code is the same as T-Lt,, i.e., the Hilbert space of 
equivalence classes of colored ribbon graphs embedded in S. 
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(a) (b) (c) 

Figure 10: The neighborhood of a vertex p in a triangulated surface can be oriented as in (a), in 
turn orienting the trivalent dual graph in (b). Attaching a "blister" above p gives (c) when r = 6. 




(a) (b) 

Figure 11: For two adjacent vertices p and p', the products Bp/Bp and BpBpi correspond to different 
internal triangulations of the same bipyramid. Hence Bp and Bp/ commute. 



7.3 Local stabilizers for the Turaev-Viro code: The Levin- Wen Hamiltonian 

Consider a vertex p, with degree r in the triangulation, as in Figure 10(a). Place a "blister" put 
on top of this configuration by introducing a new vertex q, with the new edges on the boundary 
labeled k' = {k[, . . . , A;',), as shown in Figure 10(c) (the edge from p to q is oriented upwards). The 
blister is a triangulation of a ball B with boundary edges labeled {k,k',m). Define 

Bp:= J2 ^^==TVc(S,(^,fc',m))|fc')(^|®|m)(m| , (85) 

a linear map on the space (C™)'^^'" of labelings {k,rh). Similarly as before, since gluing two balls 
together results in another ball, Bp is a projection. The additional factors compared to Eq. (84) 
are necessary because the outer edges with labels m remain on the boundary when gluing, and the 
number of vertices on the boundary stays the same. 

For different vertices p and p', the operators Bp and Bpi, acting on (C™)®'^', commute. This 
is clear if p and p' are nonadjacent. When p is adjacent to p' , both products BpiBp and BpBpi 
correspond to attaching the same bipyramid to S, only with different internal triangulations, as 
shown in Figure 11 in the case that p and p' are connected by a single edge. Since the Turaev-Viro 
invariant does not depend on the internal triangulation, we conclude BpBpi = Bp' Bp. 

Similarly, the product J^^ Bp of over all vertices of the triangulation corresponds to a "thicken- 
ing" of S, and therefore 

[]i3p = TV(Sx [-1,1]) . (86) 

p 
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Hence the Turaev-Viro code equals the simultaneous +l-eigenspace of the Bp projections. 

In fact, these Bp operators are the same as the Bp plaquette operators in the definition of the 
Levin- Wen Hamiltonian (41). To see this, evaluate the invariant TVc(-B, {k,k',rn)) associated with 
Figure 10(c) explicitly. With fixed boundary labels (k,k',m), the only free index is the label i of 
the internal edge {p,q}- A typical tetrahedron has the form shown in Figure 9(c), whence 

4, ^kt-ii*kt iK-i(K)* /o^N 

9t = I = I ■ (87) 



Substituting into Eqs. (83) and (85), we obtain 



k,k',rn * ^*=1 ^ 

Considering the dual graph, in Figure 10(b), we recognize in Eq. (88) the definition (80) of the 
Levin- Wen Hamiltonian's plaquette operators. 
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A Ribbon graphs for general categories 

In this appendix, we collect various facts about ribbon graphs and the fusion basis states discussed 
in the main text. In Appendix A.l, we summarize a few general identities satisfied by ribbon graphs. 
In Appendix A. 2, we give alternative definitions of the inner product on the ribbon graph Hilbert 
space 7i^, and use this to show that the fusion basis states are mutually orthogonal. Appendix A. 3 
provides a few explicit examples of fusion basis states for the Fibonacci model. We conclude with 
a short explanation of the generalization to higher-genus surfaces in Appendix A. 4. 

A.l Ribbon graph identities 

Here we state a few useful equivalence relations between ribbon graphs. First, bubbles can be 
removed by the rule 




Next we summarize some of the main properties of the vacuum lines defined in Eq. (71). 
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Lemma A.l. The following identities hold, regardless of the contents of the shaded region: 




' 'J ) =v 'J 



V6 



Proof. To prove Eq. (90), we apply an F-move, getting 



j 



(i 



di 



djdk 



By symmetry, this equals the right-hand side of Eq. (90). Eq. (91) also follows: 

j j j 

For the proof of Eq. (92), observe that necessarily 



(90) 
(91) 

(92) 



(93) 



(94) 



a 



(95) 



for some scalar aij E C. Taking the trace, i.e., closing the loop with label j, gives aij = DSij/dj. 
The claim then immediately follows because 



1 ^ , (±>i 



V 



where the parenthesized term is equal to 6ji by the unitarity of the ^-matrix. 
An immediate consequence of Eq. (92) in Lemma A.l is 



(96) 



□ 



(97) 
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Indeed, applying an F-move to the edge with label i gives 

i i' i 






■I k 

J 3 3 



3 f 



(98) 



I)6iif5jj>Fj 



3 



and the claim follows by Eq. (59). 

A. 2 The anyonic fusion basis states form an orthonormal basis 

Fix a pants decomposition, corresponding to a rooted binary tree T, of the n-punctured sphere 
S = S„. In this section, we show that the anyonic fusion basis states defined in Section 6.5 form a 
complete, orthonormal basis. 

The proof is straightforward, but is made inconvenient by the definition of the inner product 
in Section 6.4. Recall that an anyonic fusion basis state is specified by two edge-labelings of T by 
particles from the category C, together with particle labels for each of the n boundary conditions. 
The corresponding state is given by placing the two fusion diagrams on T, x {±1} C T, x [—1,1], 
adding vacuum loops around n — 1 of the punctures, closing up the boundary conditions, and then 
reducing to a ribbon graph in S. For example, Eq. (78) shows a typical fusion basis state on S4, 
with the pants decomposition before reducing to two dimensions. 

To compute the inner product, a ribbon graph in S x [—1, 1] must first be reduced to E, and 
then to an orthonormal basis of ^s, as for example in Eq. (70). To avoid these steps, we will 
define a new inner product, the trace inner product, that can be computed directly on ribbon 
graphs without first reducing them to two dimensions. Although the definition is quite different, 
the trace inner product {■\-)tr agrees with the original inner product (•!•) from Section 6.4. It is 
also convenient for extending the definition of "Hs to more general surfaces (see Appendix A. 4). 

Slightly more background on ribbon diagrams is necessary. A line in a two-dimensional ribbon 
diagram represents a narrow ribbon that is flat in the plane. A full twist in the ribbon can be 
drawn as a kink, as in Eq. (63), and is equivalent to adding a phase to the untwisted ribbon. The 
same is not true for half twists, though, and the ribbon-graph notation we have developed so far 
does not allow for specifying half twists. For all the equivalence relations it should be assumed that 
the involved ribbons all share the same up/down orientation. 

Let S be a compact, orientable surface, with a marked point on each boundary component. 
Consider two ribbon graphs in S. To compute their trace inner product, embed the two diagrams 
in S X {±1} C S X [—1, 1]. Conjugate all particle labels in the first diagram, and also give them the 
outward orientation, i.e., turn all ribbons upside down. Next, connect together the ribbons at the 
boundary points. (Since the ribbons on one half have been turned upside down, there is no need 
to twist the ribbons to connect them together; a twist would have introduced a phase ambiguity.) 
Finally, reduce the resulting ribbon graph to two dimensions. In order to carry out this reduction, 
it may be useful to introduce pairs of half-twists on some ribbons in order to flip their orientation 
in local regions, so that the equivalence relations can be applied. Now, the inner product is the 
coefficient of the empty or vacuum diagram, divided by the product over all boundary labels i 
of \/di. 
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Let us first sketch the argument that the trace inner product agrees witli the original definition 
of the inner product. For computational basis states on the annulus S2, we have 




Vdkdi 




(99) 



where on the right-hand side the two ribbons labeled i should be matched up and then the dia- 
gram reduced to its vacuum coefficient. Carrying out this calculation — beginning with an i^-move 
between the two ribbons j and j' — indeed gives Skk'^w^u'^jj' , agreeing with the standard inner 
product. The general calculation of the trace inner product on computational basis states on the 
n-punctured sphere S„ is quite similar; for example, 



ai a2 as a'l a'2 a'^ 



CiNT '^IX^/ ' Jda-^daida^d, 
C2 Co 




(100) 



which by Eq. (89) simplifies to S^^,5^g. 

Next we argue that the anyons are orthonormal under the trace inner product. Again we begin 
with the case of the annulus S2. The anyonic fusion basis states on S2 are the ribbon graphs 
k^e{i,j) of Eq. (77) with 5ij(,* = 6i>j'e>* = 1. 



Lemma A. 2. Under stacking, 

I'^k'ii'j') k'^eihj) = Sii'6jj>6kk' T^y 
Proof. By Eq. (90) in Lemma A.l and Eq. (97), 



^ dk 
didj 



(101) 





V6iii5jj 



(102) 



□ 



In particular, "^^k^ k{h j) is an idempotent corresponding to a realization of the fusion space V^^ 



Now the trace inner product of e'^k'i'i' if) and k'^iihj) is given by l/y/dyTi times the vacuum 
coefficient of their stack i'^k'i'i' ^j') k^tihj)^ where additionally the ribbons £ and i' are connected 
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together. The coefficient of the trivial particle 1 in a vacuum line is di/T) = 1/T), so by Lemma A. 2, 

1 




Sii'6jji6kk'Su' 



(103) 



Define anyonic fusion basis states for Tis^ 



fc^Xa, b, c, d) = «i 




(104) 



Analogous states k'^ ^a^h^c^d) G "H^^ can be defined for any n > 3 using a "standard" pants 
decomposition of S„. We show that these states, with fusion constraints satisfied, are mutually 
orthogonal and unit-normalized. 

Lemma A. 3. Let ^' = k^ ^3,b,c,d) and \&' = k'^ ^lio-' ,c' ,d') be the states in « > 3, 

as introduced in Eq. (104)- Then 

{qf'\qi) = dkk'^gi^dw^^i^cc'^dd' ■ (^0^) 
Proof. Although the calculation is fully general, we continue to illustrate the case of n = 3. Ap- 
plying F-moves, we find 




(106) 



and similarly for ^' . Thus we can separate the calculation of the inner product (^''|\&) into two 
parts, the inner products of the topmost states in '11^,2 ^"^^ the inner product of the lower doubled 
trees. From Eq. (103), this gives 

-,a' a'* 1 \ * ^a,- a* 1 



(107) 



fh^fh! i 



where 




(108) 



39 



and is defined similarly. 




Completeness of anyonic fusion basis for ^s„ 

It remains to show that the anyonic fusion basis spans "Hs^- We do so by verifying that the number 
of standard computational basis states equals the number of anyonic fusion basis. 

The proof is by induction in n. For n G {0, 1}, is one-dimensional. For n = 2, E„ the 
annulus, with boundary labels k* and i, the numbers of computational basis states and of anyonic 
fusion basis states are both 

^6ki''j*6ije* . (110) 
hj 

See Eqs. (70) and (77). 

For T-Lj:„ with a fixed root puncture, let Sn{k) be the number of standard basis elements with 
label k on the root. This function satisfies the recursion, for any m, 2 < m < n — 1, 

Snik) = '^Sijk*Sm{i)Sn+l-mij) ■ (HI) 

Similarly, let an{k) be the number of anyonic fusion states in "Hsn with root labeled k. We will 
show that an{k) satisfies the same recursion, Eq. (Ill), as Sn{k). For this purpose, let a„(^, £') be 
the number doubled fusion diagrams on a tree with n leaves, with arbitrary boundary-consistent 
labels on the n — 1 leaves and the doubled anyon label at the root. Since picking a boundary- 
consistent label k for the root puncture gives a valid anyonic fusion basis state, we have 

anik) = Y,dnii,i')5u'k* . (112) 
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Moreover, d satisfies the recursion, for 2 < m < n — 1 




(113) 



Reinserting this into Eq. (112) gives 





= ^ 5«'fc*«m(^)an+l_m(^') 

e,e',k 



(114) 



the desired recursion. Here in the second step we used 




(115) 



which is a consequence of the associativity of fusion, Eq. (56). In particular, the left- and right-hand 
sides respectively equal the number of ribbon graphs 



with fixed boundary conditions ,j,j' ,k*). In other words, both sides of Eq. (115) count the 
dimension of the space of ribbon graphs embedded in a ball, with five fixed boundary labels. 

A. 3 Application to the Fibonacci model 

Let us now specialize these statements to the Fibonacci model, and give a few examples of the 
anyonic fusion basis states described in Sections 3.1 and 6.5. Here the F and R tensors are specified 
by Eqs. (35) and (20), respectively. Crossings can therefore be resolved according to Eq. (19). 

Computation of idempotents on T^gj 

For the annulus S2, anyonic fusion basis states are equivalent to idempotents, up to normalization, 
when the boundary labels are identical. We therefore restrict our attention to the idempotents. 
According to Lemma A. 2, the idempotents can be written as 




and 




(116) 
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It = - 



V 



(117) 
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where D = Vl + t^. Writing out this definition and resolving the crossings, it can be verified that 
these expressions coincide with the idempotents given in Eq. (24). For example, we get 



p2 ir 



+ T 




+ re 



-37ri/5 




+ re 



37ri/5 




(118) 



For later reference, note that the same type of calculation gives 

Q 



-37ri/10 




\ 



^Svri/lO 




(119) 



Computation of anyonic fusion basis states for 'Hy;^ 

For the 3-punctured sphere, we have for example the following three fusion basis states. 



TT t1 

X = 

Ir 




X 



Tl 





(120) 



where the subscript i indicates that we choose to have an edge ending on the boundary. 

Let us reduce to two dimensions the first diagram. Apply an F-move to the left leg to get 



TT t1 1 

It 



^6' 



I 



+ 




\ 

\ 
I 



:i2i) 



Removing the twists with Eq. (20) and inserting from Eq. (24) the expressions for the ttq and rl 
idempotents gives 



TT t1 

X = 

It 



-47ri/5 
diTi/lO 



1 



O 







o 



+ 



(122) 



Since (9|(p) = 1/r, this state indeed has unit norm, as asserted by Lemma A. 3. 
A. 4 Generalization to higher-genus surfaces 

Our results extend to a surface S of higher genus. To describe a fusion diagram basis of S, fix 
a handle decomposition of S into an n-punctured sphere with handles attached to some of the 
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holes. A handle is a punctured torus, shown in Figure 12(a). An anyonic fusion basis for 1-Ly. is 
defined by extending a doubled anyon fusion diagram in $]„ x [—1, 1] by attaching at each handle 
one of the diagrams in Figure 12(b). Here o+a_ is the label on the leaf, and (a+,6+,6^) and 
(a_ ,h-,h*_) are fusion-consistent triples from C+ and C_ , respectively. These two choices result in 
bases that are diagonal with respect to Dehn twists along the loops {a,/3} and {0,7}, respectively, 
in Figure 12(a). 

Together with F-moves, the change of basis matrix 5"^+'^- = {S^^i^y' ^ ^ ) defined by 




(123) 



fully specifies the action of the mapping class group of S. (For the torus, there are two inequivalent 
bases, and a similar change of basis matrix S.) By similar arguments as those used in Section 6.5 
for the i?-matrix, one can derive that Sti^u~ 1.1.= -^tt' i. B?,~ , , where 

' 0,0 ,0+0- o',,b+ b ,o_ ' 



b h' b b' 




(124) 



We refer to [Kit06, Appendix E] for a proof of the unitarity of these matrices. 



B Discretizing and Gluing: Proofs 

In this appendix, we prove Lemmas 4.1 and 4.2. Lemma 4.1 states that the simultaneous +1- 
eigenspace 'H^'^ of all plaquette- and vertex-operators of the Levin- Wen Hamiltonian (41) is iso- 
morphic to the ribbon graph space H^. In preparation for the proof of this fact, we derive the 
following auxiliary statement. 

Lemma B.l. Consider two surfaces E and T,' , where T,' is the same as S but with a puncture 

(£ 1) 

inserted at p E S. Let i be a labeling of the boundary points of S, and let H-^l be the space of 
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ribbon graphs on S' with open boundary condition on p. Consider a state {"^t,) G We deform 

the ribbon graph locally around p, and regard the resulting ribbon graph as an element \ G 'H^/'^'*- 

Let B : wjf'^^ ^ be the map which adds a vacuum loop around p, divided by D. Then i?|<I>s') 
is independent of the initial deformation (i.e., l^s'))- Furthermore, the map 

A • ^ ^E' (125) 

is a norm-preserving isomorphism. 

Proof. Observe that the state i?|<I>s') has a vacuum loop inserted around the puncture p. Since 
locaUy deformed configurations of ribbon graphs near a puncture enclosed by a vacuum loop are 
equivalent according to Lemma A.l (90), the state B\^Y:r) indeed does not depend on the initially 
chosen deformation. This implies that the map A-.n^j.^ -H^V^ is well-defined. It is easy to see 
that the map is surjective; a preimage of an element |^'.^) G Ti^;^'^ is simply the same ribbon graph 
embedded in S (this can be checked with Lemma A.l (91)). 

Now consider the inner product (^e|^s) of two states |^s)>|^s) G ^|]- Since local defor- 
mations do not change the equivalence class, we may assume that these ribbon graphs avoid p. 
To evaluate the inner product (^'s'l^s') of the images \^y:') = M^s), |^s') = A|^x;); consider 
a ribbon graph basis B' of "H^l^^ whose ribbon graphs are obtained by attaching a tadpole with 
head surrounding p (and fusion-consistent labels) to the elements of a ribbon graph basis B of 
Ti^. Evaluating the inner product (^^'I^E') with respect to this basis reduces to the evaluation of 
(^fsl^E) and the norm of a vacuum line around p. This gives 

(§E'|^E') = ('I'eI^s) , (126) 

which implies the claim. □ 

The proof of Lemma 4.1 is now immediate. 

Proof of Lemma 4.1. Let |^'s) G 71^,, and let |^*a) G Ti.)^ be a state obtained by locally de- 
forming I^e) to avoid the points in V. Observe that the state B\^/s)^ where B is the product of 
all plaquette operators (cf. (42)), has vacuum loops inserted around each puncture p G V. This 
corresponds to the situation described in Lemma B.l. By inductively applying this lemma to all 
p £ V, we obtain the claim. □ 

We now prove the following more general version of Lemma 4.2, which describes the effect of 
gluing. 

Lemma B.2. Let T, Ta, Tb, b^, 6_ and e be as in Lemma 4. 2. Then the anyonic fusion basis 
state \i, d)j- is given by 



I A d)f = J2J ^^Kb-n* dA)f^ 14, dB)f^ \n)e (127) 
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Proof. Consider the off-lattice ribbon graph \l,d)j] G Tij] described in terms of a ribbon graph on 
S X [—1,1]. By doubhng a vacuum loop around a boundary component twice, and pulling out the 
loops, we can obtain two vacuum loops running parallel to the gluing curve 7 if we introduce an 
additional factor of Applying an F-move then results in 




(128) 



We can analyze the effect of applying the plaquette operator Bq in the off-lattice picture: it adds 
a puncture, a vacuum loop around it and a factor of i. It transforms this state to 




(129) 



where we used property (90) of vacuum lines. With Eq. (97), we conclude that this is equal to 




(130) 



It is easy to see that applying the product OpT^g of remaining plaquette operators results in 
the state on the right-hand side of Eq. (127). □ 

C The Witten-Reshetikhin-Turaev invariant 

In Section A. 4, we sketched how, based on a modular category C, the Turaev-Viro construction 
yields a representation pdc,y, of the mapping class group MCGe of a surface S. Furthermore, this 
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Figure 13: The anyonic basis state |fc,g) £ T~ic,dMg, here for g = 4. All edges carry the trivial 
label 1 of C. 

representation is described by the doubled category DC = C <^ C* . More generally, any modular 
category C gives rise to a (projective) representation 

PCS : MCGs ^ GL(?^c,s)/ < e2-^^/24 > ^ ^^3^) 

where Hce is a Hilbert space of anyonic fusion diagrams. (In this expression, < g*^'^*'^/^'^ > denotes 
the cyclic group generated by g^'^*'^/^^'^!, where c is a scalar called central charge.) A certain 
matrix element of the representation (131) for the boundary S = dMg of the g-handlebody Mg 
(cf. Figure 14) defines the Witten-Reshetikhin-Turaev (WRT) invariant [Wit89, RT91], as we now 
explain. 

A closed, oriented 3-manifold M can be represented by a Heegaard splitting {g,x), i.e., a 
genus g G N and an element x of the mapping class group MCG{dMg) of the surface dMg. Roughly, 
the element x specifies how to glue two copies of Mg together to obtain M. The pair {g, x) 
uniquely specifies the equivalence class of M under homeomorphisms, up to (i) a certain stabilization 
move {g, x) — t- {g+1, x) which corresponds to attaching a 3-sphere with a standard genus-1 Heegaard 
splitting, and (ii) multiplication of x by the subgroup MCGg^^^^ of mapping class group elements 
which extend to homeomorphisms of Mg. 

For a modular category C, the Witten-Reshetikhin-Turaev invariant of a manifold M with 
Heegaard splitting {g, x) is given by 

\NRJc{M):=V9-^{vc,g\pc,9{x)\vc,g) , (132) 

where |wc,g) is the unit-normalized vector corresponding to the anyon diagram shown in Figure 13. 
Invariance follows from the fact that the vector |fc,g) is invariant under the action of MCGqj^^, and 
the factor D^^i compensates for stabilization moves. Eq. (132) is the Crane-Kohno-Kontsevich [Kon88, 
Cra91, Koh92] version of the WRT invariant. Its original [RT91] and more commonly known defi- 
nition is based on the Dehn surgery presentation of manifolds. Piunikhin [Piu93] (see also [Koh02, 
Section 2.4]) showed that these definitions are equivalent (for C = SU{2)k)- 

For a closed 3-manifold M and a modular category C, definition (83) of the Turaev-Viro invariant 
reduces to 

tVc(m) = p-2|v'mI J2 U^He) n at- (133) 

labelings (p ^ tctrahedra t 

The invariants (132) and (133) are related by the categorical double, that is, 

TVc(Af) = WRTdc(M) . (134) 

Eq. (134) was shown by Walker and Turaev [Wal91, Tur94] (see also [Rob95]) and holds more 
generally if C is a spherical category [TurlO]. As explained in Section 6, the double DC takes the 
form C i^iC* for a modular category C and (134) becomes 

TVc(M) = |WRTc(M)p . (135) 
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This is because pc(^c* — Pc® PC* is a tensor-product representation of C and its conjugate, the state 
of Figure 13 factorizes as \vc(^c*,g) — l^c,g) ® \vc*,g), and the total quantum dimension of C ® C* 
is equal to V^. Eq. (135) has direct application to quantum computing: it facilitates the proof 
in [AJKRIO] that approximating the Turaev-Viro invariant is a BQP-complete problem. 

In the remainder of this section, we sketch a proof of (134) (for modular C) based on our 
description of doubled anyonic fusion basis states arising in the Turaev-Viro code. 

A first step is the observation that, by definition of a Heegaard splitting, a manifold described 
by {g,x) is the result of gluing together three manifolds: two copies Mg, M'g of the (7-handlebody, 
and a "mapping cylinder" or cobordism defined by the mapping class group element x. Contracting 
the Turaev-Viro-tensor networks corresponding to the handlebodies gives two states \Z{Mg)) and 
\Z{M'g)) on T-LdMg-, and T-idM'g-, respectively, while the mapping cylinder gives rise to a linear map 
PC^C-ix) ■■ TidMg T-LdM'g- The state \Z{Mg)) is given by 

\Z{Mg)):=^JVc{Mg,x)\x) (136) 

X 

where the sum is over labelings x of the edges of the triangulated surface dMg. (In the context 
of TQFTs, such states associated to 3-manifolds are sometimes called partition functions [Wal91].) 
Because of its definition in terms of the Turaev-Viro invariant, this vector is in the Turaev-Viro 
code subspace on dMg, defined by the projection (84). A similar expression holds for \Z{Mg)). 
The linear map pc^c*ix) can be interpreted as being composed of F- moves (i.e., local changes of 
the triangulation) because of its definition in terms of contractions of the F-tensor, and pc^c* is a 
representation of MCG{dMg). This is the representation (131) arising from the doubled category 
C* (cf. Appendix A.4). 

Combining these three components, we conclude that the Turaev-Viro invariant has the form 

JVciM) = {Z{M'g)\pc^c*{x)\Z{Mg)) , (137) 

which already bears some resemblance with the WRT invariant (132). Indeed, Eq. (134) follows 
from (137) and the following lemma, which determines the anyonic fusion basis state corresponding 

to \Z{Mg)). 

Lemma C.l. Consider a genus-g handlebody Mg with triangulated boundary dMg, and let \Z{Mg)) 
he the Turaev-Viro codeword given by Eq. (136). 

(i) The off-lattice ribbon graph underlying \Z{Mg)) (cf. Lemma 4-1) consists of a vacuum loop 
around each handle, as in Figure 14- 

(ii) The state \Z{Mg)) has squared norm \\\Z{Mg))f = 

In particular, up to a phase, \Z{Mg)) is equal to 'D'-^^'^\vcr^c* ,g) j where |uc®C*,s) is the normalized 
anyonic fusion basis state given in Figure 13. 

The proof of this lemma proceeds by induction in the genus using the fact that a genus-^ 
handlebody can be built from a genus- (5 — 1) handlebody by gluing together two discs from the 
boundary, as illustrated in Figure 15. To derive a gluing formula describing the transition from 
\Z{Mg^i)) to \Z{Mg)), we may locally choose a convenient triangulation, and consider the result of 
contracting the corresponding tensor indices against each other. For example, we may choose each 
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Figure 14: The state \Z{Mg)), here for (7 = 4, is the projection of this ribbon graph, with a vacuum 
loop around each handle. 




Figure 15: Two surfaces can be joined by puncturing each and identifying the boundaries of the 
holes, as in Section 5.1. In the case that the surfaces are the same, the boundary of a 3-manifold, 
gluing together two disks from the boundary adds a handle to the manifold. The Turaev-Viro state 
for the resulting manifold has a vacuum loop going around the handle. 

disk to be a single triangle, and moreover may take the triangle to be degenerate, with two of its 
sides the same edge. In this case, the dual graph is locally a tadpole. From the characterization of 
the Levin- Wen plaquette operators in Section 6.6 each as adding a vacuum loop around a puncture 
within the plaquette, and Eqs. (90) and (68), the qudit along the tail of a tadpole is fixed to 1. 
Hence the state \Z{Mg-i)) can be factored as 



where the second term corresponds to the two tadpole tails, the third term is the state of the 
tadpole heads, and \v') is a state on the remaining qudits. This expression then implies the desired 
result for \Z{Mg)) when using the definition of TV^. Below we will give an alternative proof that 
avoids using degenerate triangulations. 

Proof. We first prove (i) by induction in the genus g. For the genus-zero case, the Turaev-Viro 
codespace on the sphere is one-dimensional, and given by the projection of the "empty" ribbon 
graph; this is the claim (i). For the induction step, we use the gluing procedure illustrated in 
Figure 15. This is equivalent to taking two patches of surface with the free tensor indices all 1, 
adding a layer of tetrahedra to project to the Turaev-Viro code, and then identifying the indices 
for patches on either side. Therefore, the new tensor that we are contracting on corresponds to a 
sphere with two separated disks of indices all fixed to 1 . By retriangulating, we may take this sphere 
to consist of just three tetrahedra, as in Figure 16(a). We want to understand this tensor with the 
six indices marked * fixed to 1. This leaves six free indices, marked . By the fusion constraints 
on the associated ribbon graph, these indices must all match. It remains only to determine the 
coefficient of this shared index; call it i. 




(138) 
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Up to a power of D, we can read off from Eqs. (83) and (81) the coefficient as 

6 F^^^ F^^^ 

where the first term, \/df , comes from the term HeeSM \Jd<f{e) in the Turaev-Viro expansion, and 



the other terms are the tensors for the three tetrahedra. Now substitute in -Fj^*^ = \j'd^.^ijk-, so 

FI^\ = 1 and Fll*l = 1/di. We find that the coefficient is proportional to di, as claimed. 

For the proof of (ii), note that the squared norm of the vector \Z{M)) associated with a 3- 
manifold M is equal to the Turaev-Viro invariant of the closed 3-manifold obtained by gluing M to 
itself. In particular, the squared norm is independent of the triangulation of M. We can therefore 
again choose convenient triangulations. 

As the ball is triangulated by a single tetrahedron, as in Eq. (81), we immediately obtain 



|||Z(Mo))f = TVc(Mo)=p-« ^ \F2ydidjdkd, . (140) 

ijklmn 

By unitarity, Y.n \Ke*Z\'^ = ^i'jm*5kmt*, so applying Eq. (57) repeatedly gives 



\\\mo))\? = ^ . (141) 

for the genus-0 case. 

To compute the norm of \Z[Mg)) for g > 1, fix a triangulation of Mg_i and consider the 
triangulation of Mg obtained by connecting two triangular faces on dMg-i using two tetrahedra 
as shown in Figure 16(b). In this figure, the face with edge labels {m,j,i) is glued to a triangle on 
dMg-i, while the face with edge labels {i, q, s) is glued to another such triangle. We claim that 

\\\Z{Mg))f=VmZ{Mg.,))f . (142) 

Using (141), the claim (ii) immediately follows from (142). 

To show (142), observe that the contribution to TV(M, x) of a labeling x locally matching 
Figure 16(b) is given by 

Tpi*jm rpq" 



kl*n nl*s 
V dmdn V dids 



yjdidj dk de dm dn dq dr ds , (143) 
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since all edges are on dM. On the other hand, the restriction x' = x|A/g_i gives a contribution 

dmdidj ■ dqd^ds (144) 

to TV(x',Mg_i). We can therefore write 



\ZiMg)) =VW\ZiMg.i)) where W= ^ ^^e'lKl* nl^ jkemnqrs)i{ijm\ (g) {qis\) 

ijkimnqrs mis 

(145) 

In (145), the additional factor of D accounts for the fact that the number of vertices on dM is 
equal to \VgM\ = \VgM'\ + |VaM"| — 1 since the two faces share a vertex after gluing. Using the 
tetrahedral symmetry of the i^-tensor, we get 

= E \T.dkKs:^\' \Fllrn i\ijm){ijm\ \q£s){qis\) . (146) 

tjmgls \knr / 

' V ' 

=:X(ijmq£s) 

We claim that 

X{ijmqis) = 5^*s*q^i*jm* for all (ijmqis) , (147) 

that is, W is an isometry on the span of valid labelings of dMg^i. This proves the claim (142) and 
concludes the proof of (ii). Indeed, we have, again by unitarity of the basis change, 

\Fqs*r = ^ni*eSe*s*q , (148) 

r 

and therefore 

\Fqs*r P l-^fcte™!^ = <J£*s*g'^i*jm*'5mM ■ (149) 

nr 

Multiplying by -j^^ and summing over k using (57) gives (147). 

Having shown (i) and (ii), we can reexpress \Z(Mg)) in terms of the anyonic fusion basis as 
follows. First, we use properties (90) and (91) of Lemma A.l to add vacuum lines around the g — 1 
constrictions between the handles in Figure 14. Since the resulting vacuum lines correspond to a 
pants decomposition of dMg, we can use the correspondence between doubled fusion basis states 
and ribbon graphs to conclude that the state \Z(Mg)) is proportional to the anyonic fusion basis 
state \vc^c*,g)- The proportionality constant is determined by the normalization (ii). □ 
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